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We propose a new and effective method to find plasma oscillatory and wave modes. It implies 
searching a pair of poles of two-dimensional (in coordinate x and time t) Laplace transform of self- 
consistent plasma electric field E(x,t) — > E P1P2 , where pi = —iu), p2 = ik are Laplace transform 
parameters, that is determining a pair of zeros of the following equation 



E P1P2 

This kind of conditional equation for searching double poles of E P1P2 (which correspond to terms 
of the type -, , M- j t^jy in double Laurent expansion of E P1P2 ) we call "general dispersion 

(Pl-Pj" j [P2-P 2 ) 

equation", so far as it is used to find the pair values {(J^ n \ k^), n = 1, 2, . . .. It differs basically 
from the classic dispersion equation ei(uj,k) =0 (and is not its generalization), where q is longi- 
tudinal dielectric susceptibility, its analytical formula being derived according to Landau analytical 
continuation. In distinction to ei, which is completely plasma characteristic, the function E P1P2 is 
defined by initial and boundary conditions and allows one to find all the variety of asymptotical 
plasma modes for each concrete plasma problem. In this paper we demonstrate some possibilities 
of applying this method to the simplest cases of collisionless ion-electron plasma and to electron 
plasma with collisions described by a collision-relaxation term —vf^- 1 '. 

PACS numbers: 52.25 Dg; 52.35 Fp. 
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I. INTRODUCTION 



Up to date the textbook approach to determining dispersion dependence of oscillation frequency on wave number 
uj(k) implies solving some dispersion equation, the latter being found using Landau rule of bypassing poles in his 
theory of collisionless damping ( "Landau damping" , seea) . This theory is based on analytical continuation of complex 
function defined as Cauchy integral 

oo 

= r 

J Z-LU 

— oo 

where the integration is performed primarily along the real z-axis. This integral is well defined for Im (uj) > and 
Im (uj) < 0. Landau considers as physically defined just the function F + (uj) in the upper half-plane of uj due to 
adiabatic in time switching on. Then he tries to determine analytical continuation of this function to the region 
Im (uj) < through the point Im (uj) = on the real axis, with deformation of the integration contour to the lower 
half-plane. 

One could say that analytical continuation through integration contour A (i.e. the real axis R) does not exist by 
its very definition as a regular function, since Ff(u) = F + (uj) is indefinite and not regular for uj £ A, i.e. for uj on 
the real axis, discriminated by physics of the problem. Suppose the fixed integration contour A' is parallel to and 
lies above the contour A. When uj lies above A' according to Cauchy theorem we have (uj) — Fa>{uj). But Fa>(uj) 
is indefinite for uj £ A' and therefore both functions F^ and Fa> for the first sight cannot be regarded as analytical 
continuation of one to another. 

Nevertheless one can resolve above-mentioned paradox and construct the analytical continuation by the following 
way. Let contour A' to lie over A by ib, b > 0. According to Cauchy theorem for Im (uj) > b we have Fa> (oj) = F^(uj), 
so 
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4>(z)dz f 4>{z + ib)dz 

z + ib-u ' [) 

where 4>(z) is supposed to be analytical in the stripe < Im (z) < b. Analytical continuation of Fa>(lo) into the region 
< Im (uj) < b can be obtained by Eq.(||) with some arbitrary b < Im (a;). In the same way if Im (uj) < one should 
use Eq.(|2) with some arbitrary b < Im (uj). In all cases we obtain the single analytic function (due to singleness 
theorem). H This function is obtained by deforming integration contour in such a way that point uj lies always above 
the contour (this is due to the + -type of considered analytical function). For Maxwellian plasma <f>{z) ~ zexp (-az 2 ) 
with |^(z)| — > oo at Im(z) — > ±oo and it has no poles in the complex plane z G C. Therefore the corresponding 
functions F (uj) are analytic functions in the whole complex plane. 

Using analytical continuation of F + (uj) in the lower half-plane (instead of F~(uj)) is closely related with the 
absence of any solutions of dispersion equation for any S in to = ujq ± iS, luq ^ 0, when dispersion equation is derived 
by substitution of plane wave solution exp(— iujt + ikx) into Vlasov equations for electrons (that is linearized kinetic 
equation + Poisson equation) in the case of background Maxwellian plasmata. The theory of Landau damping can also 
be considereda as an equivalent to solving the analytically continued to the region Im (w) < dispersion equation, 
that is equation F + (oj, k) = 1 with analytically continued left-hand-side function F + (uj, k). 

The mathematical part of Landau theory (see f.e.LJ) contrary to the said inOu does not give rise to any objections. 
However there are still the following logical objections: 

1. Landau solution in the form of a damping wave must satisfy primary Vlasov equations, but it does not. 

2. In the original Landau contour integral F(uj) there is the physically discriminated real axis of integration over 
velocity component v x . One solves some initial problem in the time interval < t < oo, and coordinate 
dependence is introduced through the common factor exp(ikx) for all times t. Here initial conditions really are 
in no connection with unknown conditions at t < 0. Writing solution in the form of inverse Laplace transform 
("original") with parameter p — —iuj physically leads one to search in the first place some finite (damping or 
oscillatory non-damping) solutions with u> = ujq — iS, ujq > and 5 > 0. These solutions correspond to the lower 
half-plane of complex frequency uj. 

The function F~(uo) is analytical in the lower half-plane, and there is no logical necessity to use analytical 
continued from the upper half-plane function F + (lo), or consider some analytical continuation into upper half- 
plane, where oscillatory solutions are exponentially divergent. As it is also well-known, the integrand of inverse 
Laplace transform L(F) in this case has no poles and, correspondingly, no solutions in the form of damping 
wave (that means the absence of solutions of equivalent dispersion equation). 



3. The application of two-dimensional Laplace transformationaa (see also this paper) shows that the asymptotic 
solution of Vlasov equations with Maxwell distribution function exists only in the form of a system of coupled 



1 For the analytical continuation of F + (lo) from the upper half-plane (this is marked by 'plus'-sign in contrast to analytical 
continuation of F^ (uj) from the lower half-plane, which is marked by 'minus'-sign) to the lower one we have 

oo 

F+(lu) = PV [ + i7r 0( w ) ( a t Im ( w ) = 0); 

J Z - LU 
— OO 

oo 

F + {uj) = PV [ + 2nT<t>(iu) (at Im (uj) < 0); 

J Z -LU 

— OO 

oo 

F-(u)=W [ (atlm(w)<0); 
J z ~u 

— OO 

F + (uj) = F' (uj) + 2in(j)(uj) (for all uj). 

Here F + (uj) is regular atiiall uj with — bo < Im(a>) < 0, where bo is defined by the analiticity properties of 4>( z )- Cf. also the 
analogous expressions inuEl. These relations and all features can easily be demonstrated with a simple example of the function 
4>(z) = 1/ (z 2 + uj 2 ) (noted by A. P. Bakulev) 
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oscillatory modes. The absence of solutions in the form of a single wave exp ^— iiot — ikx^j , where to is complex 
value, can be proved by reductio ad absurdum: the substitution of this solution into Vlasov equations leads 
to a dispersion equation which has no solutions of the type lu = f (jcj . Landau theory for the solution in 

a more general form f(t) ■ exp (j,kxj ought to be considered as a part of the classical proof by reductio ad 

absurdum: using one-dimensional Laplace transformation for the half-plane of damping waves Im (lu) < one 
obtains, contrary to initial supposition, that damping solutions in the form f(t) ~ exp (— itut) are absent, and 
moreover the analytical continuation into the half-plane Im (to) > leads, also contrary to initial supposition, 
to exponentially growing solutions. Besides that the latters satisfy not the Vlasov equations but some other 
equations with another dispersion relation, that proves in its turn the failure of the initial supposition about 
the assumed form of solution. 

Therefore there is discrepancy between mathematical correctness of Landau theory and unjustifyness and unnecessity 
of its real application. _ 

As it has been pointed out ina the solution of "Landau problem" is really non-damping standing oscillations (but 
not a travelling damping wave) as it has been proved iroo, and the desired damping solution can be attained through 
the transition to a non-Maxwellian electron velocity distribution function ( "cut-off" velocity distribution) . Probably, 
it is equivalent to solving dispersion equation which can be solved with respect to 5 owing to using just non-Maxwellian 
background distribution function. But evidently, such solution is not completely defined due to the arbitrariness of 
Maxwell distribution "cut-off". And besides that, under special condition of "cut-off" v < c where c is the light 
velocity, one should in any case solve the relativistic kinetic equation. 

A quite different, but general way to consider plasma oscillations in frames of linearized Vlasov equations was 



proposed inEffl. It relies on two-dimensional (in coordinate x and time t) Laplace transforms f P l P2 , E P1P2 for the 
perturbation f^(x,t,v) of distribution function and the field E(x,t), where pi, P2 - are the parameters p\ = —ilu, 
P2 = ik of Laplace transformation. Depending on initial and boundary conditions one obtains different analytical 
expressions for the function E P1P2 . The equation for double poles (in p\ and p^) of this function defines different 
asymptotic oscillatory modes as the pairs (uj^ n \k^), n — 1,2, ... . So, it is natural to consider the equation for poles 
of E plP2 as "general dispersion equation" . The general asymptotic solution is obtained as a sum of exponential modes 
with coefficients defined by residues, but substitution of this solution /W (x,t,v) and E(x,t) into Vlasov equations 
does not lead, in general, to any dispersion equation connecting lu and k. That is, the dispersion equation, in common 
sense, does not exist. Indeed, it ought to suppose that, vice versa, for a given concrete problem its asymptotical 
solution in the form of coupled oscillatory modes must identically satisfy Eqns.(^)-(^), and this can be tested by 
direct substitution of solution in Eqns.(|3|)-(U). 

This approach appears to be more general, allowing one to find additional oscillatory modes also in the case of 
non-Maxwellian distribution functions, whea,the usual dispersion equation could take place and have solutions. 

These statements which were developed ina appear to form a new and unexpected, though very simple, approach to 
the problem of plasma oscillations. In this paper we demonstrate the very possibility of applying the two-dimensional 
Laplace transformation for finding ion-electron oscillations and waves in collisionless Maxwellian plasma and electron 
oscillations in low-collisional Maxwellian plasma in the case of electron-neutral collisions described by the collision- 
relaxation term —vf^\ where v ~ const is collision frequency. (But, if Landau theory is wrong, then it is necessary 
also to revise collision corrections for accurate taking into account the Coulomb collisions including one obtained in 
□ with the method of expansion into asymptotically divergent series over small parameter ^.) 

To demonstrate our approach we consider here the simplest cases of one-dimensional (in x) plane longitudinal 
plasma waves based on kinetic equations for electrons, correspondingly, ions 



d/ (le) | v d/ (le) \e\E x dfS 

dt x dx m dv x 

df^ | y dfM | \e\E x dfi 

dt x dx m dv x 



= (3) 
= (4) 



and Poisson equation 

dE x 



dx 



= -A7r\e\n e / f^d?v + 4nz\e\m / f^d 3 v; n e ~ n z . (5) 
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II. THE BOUNDARY PROBLEM OF EXCITING ION-ELECTRON OSCILLATIONS AND WAVES 



The plasma is assumed Maxwellian, homogeneous and infinite in space and time. The boundary condition is 
assumed to be given in the plane x — (the plane geometry). The one-sided Laplace transformation allows one 
to obtain absolutely definitive solution in the region < x < oo. But for the same boundary condition one can 
also obtain one-sided solution with the help of analogous Laplace transformation in the region — oo < x < (after 
substitution x' — —x). So one obtains the united solution in the whole interval — oo < x < oo. The solution will be 
continuous at x = 0. The derivative df^ jdx at x = can be either continuous or discontinuous. The same is true 
also for the one-sided Laplace time transformation at — oo < t < oo for some given initial condition at t = 0. The 
united solution will not be specially analyzed here. 

Applying two-dimensional Laplace transformation to perturbations of the electron and ion distribution functions 
/ (le) , / (H) 



l 



(2ni) 

dfW{x,t,v) 1 



■>PlP2 



dp 1 dp 2 
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(72 +i00 
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2m 



t/ (1) (0,t,v) 



e P2X d P2 , 



(G) 
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and analogous expressions for 
conditions 



d} w (x,t,v) 

at 



dE(x,t) 



E(0,t) =E e- l(it 
fW(0,t,J)=ae- i P t 



in eqns.(|^)-(^); assuming for instance oscillatory boundary 
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one obtains for electrons/ions: 



Ale) 
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Pi + V x p 2 



v x e 



kT P \ 2nkT P 
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(li) = 

P1P2 Pi + V xP2 



v x \e\ f M \ 3 ^ 2 _ju 



kT \2irkTJ 



e 2kT E. 



P1P2 



pi + ip p2 



where Maxwell functions have been used: 



( m 
\2m\T e 



3/2 



M 
2nkT 



3/2 



(10) 
(11) 

(12) 



and fpl^ (v), fpV(v) are corresponding one- dimensional (either in t or x) Laplace images of boundary and initial 
conditions (see (fil)). Then equation for finding E pip2 can be written in the following form: 



E, 



ini'2 



P2 



4ire 2 n e ( m 
kT e I 2-rtkTe 



3/2 



'2™r c . v d 3 v Aire* 



n e ( M \ 3/2 7 
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En 



47T e n, 



Pi + v x p 2 kT 

(a e - ai) (v) v x d 3 v 



Pi + i(3 

We can transform integrals in the l.h.s. of eq.(|13| 



{pi + v x p 2 ) (pi + if3) 



e ^ft ■ v x a v 
Pi + v x p 2 

—00 

00 ^2 1C) - til^) d 3 V 



(Pi + V x p2) 



(13) 
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•J x d 3 v „ [e~~--vld z v /27rfcT e V 



P2 



2 P2 / — 5V- - - 1 V 2 % 2 ' ( 14 ) 

Pi + v x p 2 J pf - Kp z 2 \ m J 2 pi - v%pl 

— oo 

where v 2 is some characteristic value of squared velocity defined by Maxwellian exponent and 7 ~ 1 (see also a 
discussion on principal value sense of such integrals incl). Note that the points pi = ^v x pi here are not really poles! 
The infinities are cancelled after accoup.t,for the analogous terms with pi ± v x p 2 in denominator from fpl (v), with 
latter being proportional to E P1P2 (seeou). Let us for simplicity omit off the consideration for a time the additive 
terms proportional to {0,t,v), (0,t,v), f (le ^ (x,0,v) and (x,0,v) and obtain 

Eg 



p _ Pa(pi+i/3) ,-.r\ 

- C/ PlP2 — , ,2 ,,2 1 K 10 ) 



p\-ap\ pf-bp: 



where uj e , uji are correspondingly Langmuir electron and ion frequencies; y/a = y/2^(kT e /m, \Jb = y/2jkT '/M are 
nearby thermal mean velocities of one-dimensional thermal movement of electrons and ions. ^ The value 7 can be 
approached in concrete calculations by means of an iteration process with more and more precise values of the poles 
(pi,P2) of expression (15[) with their subsequent substitution in eqn.(14). The equation for finding a pair of poles of 



the function E plP2 ([15 ) we call "general dispersion equation" for given boundary and initial problem. 

If the terms a e (v), a, (v) and f P l e \ fp\ are non-zero, they can be presented as sums of symmetrical and antisym- 
metrical (with respect to v x — > —v x } parts and be reduced by analogous way to integrals over interval < v x < 00 with 
the integrand denominator (p\ — v^p^)- This might lead also to an appearance of new plasma modes. The plasma 
boundary oscillations apparently can be realized by supplying periodic electrical potential to wide flat electrode grid 
immersed into homogeneous plasma volume. 

Asymptotical both in x and t plasma modes (uj^ n \ fc(™)) are defined by a pair of residues (in p\ and p%) at poles of 
function (|l5|). The case P2 = (that is, relative displacements of electron and ion parts as a whole) leads to the pole 
Pi = ii-v/wg + wf. It corresponds to the eigen-mode of non-damping plasma oscillations with the frequency (but not 
residue) independent on (3. To excite this mode the boundary field should be increased fast enough up to the value 

Eq. 

Besides that there is a mode p 2 = 0, p\ = —i(3, which corresponds to forced plasma oscillations of electrons and 
ions as a whole with the frequency of excitation source. 

The more interesting case is p\ = — i[3 and p 2 being the root of equation 

oj 2 us 2 

1_ f3 2 +a P 2 2 ~ (3 2 +bpl "° • (16) 
Its solution can be found analytically in general form: 



We do not consider here all particular cases, but instead note that at common conditions M 3> m, T e 3> Tj with 
(3 ^> LOi one obtains the following roots: 



^--^-L; p^>~±i0 X / ±1^, (18) 



so the wave speeds (if there are no standing waves) are 




-"VP 2 V M y/l-Jt/P 

that is the modes are the Langmuir ion waves. 



It is necessary to correct iJi the term £o/pi by the term -Eo/P2- In Eq.(28) the second pair of poles must be 
±i^2Ae-^A^,p 2 = 0) 
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If (3 <C u)i, the expression under the root sign can be expanded in the small parameter 



Aabf3 2 (f3 2 - uj 2 - uJf) 
At [a (/3 2 — cj? ) + b (/3 2 — cj 2 )] < one of solutions will take the form 



2 „ ~P 2 "l . (3,4) _ »/3^e , . 

P2 - „. .2 , I., .2 ' ^2 - ± 7= / 2/ J ) <^ U J 



so the travel speeds are 



^^^(M.^(.4), (21) 



that is these solutions correspond to the modes of non-damping ion-acoustic waves. The other pair of values P2 are 



ab 



and correspond to exponential damping and exponential growing modes. 

Besides that in all the cases the damping solutions are not of the specific Landau damping type. 

As it has been noted inta, the presence of growing solutions is connected with the inconsistency of the initial and 
boundary conditions on the field E(0,t) and functions /W (0, t, v) and /W (x, 0,v) , which are not independent. 
The initial and boundary conditions must be given in such a way, that the residues sums in the coefficients of 
exponentially growing modes be cancelled. Such consistency condition leads to linear equations connecting pre- 
cxponential coefficients of different modes and makes the existence of plasma oscillations not independent on electric 
field impossible. 

Laplace images of the form t^jfy or (|2^) are general solutions of given concrete problems after applying inverse 
Laplace transformation. It stays now but an open question, whether one can with integration over p2, say, to pass to 
1-dimensional Landau initial problem and to his 1-dimensional Laplace image. This is in accord with above statement 
that any plasma oscillatory problem is not strongly only initial or only boundary, but instead mixed problem with 
coupled oscillatory modes. 

It should be noted that a single-valued solution of Vlasov first-order partial differential equations for the self- 
consistent plasma electric field E con (x,t) is completely defined by initial condition E con (x,0) and by a single-point 
boundary condition E con (0,t) only. It appears to contradict to a possibility of an arbitrary setting electric field in 
several points using some electrode system. This contradiction is removed by taking into account two different parts of 
the total electric field E(x, t) = E con (x, t) + E ext (x, t) in eqns.(^-|^), where E ext (x, t) is inconsistent external field (its 
inconsistency means that it is defined by external source, not by plasma). In practice one usually transits from field to 
potential according to E ext = —d<p ex t/dx, and E con must satisfy Poisson equation (||) and is defined as self-consistent 
field with boundary condition _E" co „(0,t). 



III. PROPER ELECTRON WAVES IN THE LOW-COLLISION PLASMA 

Adding into the r.h.s. of eqn.(||) the collision term — v (v) in close analogy with the previous section we obtain 

p _ P2(pi+i/3) 

I i 7 ; — 2 

(Pi+u) -ap% 

where v (v) is an effective collision frequency between electrons and neutral particles. At p2 = one obtains the 
pole pi — ±ico e — v - that means the damping in time Langmuir oscillations. Besides that there are non-damping 
oscillations with the boundary field frequency (p2 = Q, pi = —ifi)- 

At pi — —i(3 one obtains also according to ( |23|) the equation for determining p2- 
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pl = ~[(v 2 -f3 2 +ojl)-2i(3v] . (24) 
a 

Omitting the elementary algebraic procedures of extracting the root square of a complex variable, we give the final 
result: 



/'-> = — — T -j= (VF 2 + 4/?^ 2 - Fj , (25) 



(yF 2 +4/? 2 !/ 2 -F 
where 

F ee i. 2 - f3 2 + luI 



At F < and small v <C \] ft 1 — oj 2 , f3 > uj e one obtains 



V^ 2 - ^ e 2 

At F > it is more convenient to rewrite eqn.(^5|) in the form 

1/2 



P2- — ^===T^V/3 ■ (26) 



(VF 2 + 4/3 2 ^ 2 + F) /3i/^/| 
and at v <C -\/^ 2 — /3 2 , ui e > (3 one obtains 



P2 ^ — T V /F2 J 2 "^' (27) 

\/F 2 + 4/3 2 is 2 + F) 



* K± ffi! T <£_» . (28) 

Va Vw 2 - /3 2 

In both cases one obtains the exponential growing solutions. But substitution of any exponential growing expression 
of the fosm exp(iwi + ikx) with complex values uj, k in eqns.(|^)-(^) leads to a dispersion equation which has no 
solutions! It means that at the same extent as taking into account the boundary condition F(0, t) one must also 
for the sake of consistency with E(0,t) take into account initial and boundary values of-,functions f^ 1 ^ (x, t, v), 
(x, t, v), in such a way that growing terms be cancelled (as it has been already noted ina). 

The possibility of such cancellation is provided with the same poles of growing solutions for E P1P2 and Laplace 
images of initial and boundary values of f^ 1 ^ and f^ ll K 

IV. CONCLUSIONS 

We have considered Vlasov differential equations for collisionless plasma and successfully solved them asymptotically 
by the Laplace transform method of operational calculus. 

We have proposed absolutely new, very simple and effective way of finding plasma oscillation modes: they are defined 
by the pairs (w (n) , fc (n) ), n = 1,2,..., which are determined as pairs of roots ( "double-zeros" ) of the "generalized 
dispersion equation" 

' = , (29) 



E PlP2 



where E pip2 is two-dimensional (in x and t) Laplace image of self-consistent plasma electric field E(x,t). Some 
additional plasma density oscillation modes appear to be determined from 

' ' 3, (30) 



Ale) " ' 

JP1P2 JPlPl 

but really this has no place. According to the required mutual consistency of initial and boundary values of 
ft 1 ^ (x, t, v), f 1 * 1 ^ (x, t, v) with boundary field F(0, t) (the solutions finiteness condition) these additional oscilla- 
tory modes must be connected with the electric field, thus the coefficients in the total sum of the modes have to be 
proportional to Eq. 



7 



The eqn.(^9|) is in principle different from equation commonly used in literature (see, for instance, the classic 
text-book exposition irfl'S) 

e,(w,fc) = 0, (31) 

where (according to the considered in this paper cases) ej is longitudinal dielectric susceptibility defined by only the 
intrinsic parameters of a homogeneous infinite plasma. Contrary to this approach the function E P1P2 accounts for 
initial and boundary conditions concretely. The finding E P1P2 , contrary to £/, does not require using Landau theory 
with its doubtful foundations (using analytical continuation, expansions in asymptotically divergent series (seeErQ), 
problems with the principle of causality and so on). At the same time the equation ( |29| ) allows one, contrary to 
(|3l|), to find all the oscillation modes in a given concrete problem. It reveals immediately the tight interconnection 
of oscillatory modes with the concrete conditions, that is with methods of excitation, excitation frequency, functional 
forms of initial and boundary conditions, whereas eqn. (|3l|) connects plasma oscillations and waves only with proper 
intrinsic parameters of plasma. 

Moreover, the required connection of initial and boundary conditions means that it is impossible completely to 
discriminate between a pure boundary or pure initial problems. So, to find the finite solution of a "boundary" problem 
it is certainly required to supply also non-zero initial values of distribution functions under consideration. 

It should be emphasized once again that in general case using the theory of plasma oscillations, based on suscep- 
tibility e and plasma modes equation e = 0, either is erroneous (because it is based on the dispersion equation which 
does not exist) or has rather limited range of applicability. 

In this paper we have considered only the simplest examples of application of suggested method. Besides of diverse 
generalizations one of the main problems here stays the account for Coulomb collisions, which requires cumbersome 
mathematical computations with revision of results described in literature (see, for instance,^) and based mainly on 
expansions in asymptotically divergent series in small parameter 5/loq. 
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